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1. Introduction 



connected regions with different dimensions. Typical CAD methods compute a CAD 
(represented by sample points) of R" which has a cylindrical structure and guarantees 
that the intersection of the CAD (the set of sample points) and every region defined by / 
is not empty. It is well known that the scale (the number of sample points) of CAD is 
usually much bigger than the number of regions defined by the polynomial. It is also 
the case that the scale of open CAD (see Definition 14) is usually much larger than the 
number of open regions defined by / ^ 0. 

In this paper, we introduce the concept of open weak CAD. Open weak CADs might 
not have cylindrical structure in the sense of classical CAD. However, open weak CADs 
preserve some geometric information, i.e ., they are open weak delineable (see Definition 
1). As the result, every open CAD is an open weak CAD. On the contrary, an open weak 
CAD is not necessarily an open CAD. 

An algorithm for computing open weak CADs is proposed. The key idea is to compute 
the intersection of projection factor sets produced by different projection orders. The 
resulting open weak CAD often has smaller number of sample points than open CADs. 

The algorithm can be used for computing sample points for all open connected compo¬ 
nents of / ^ 0 for a given polynomial /. It can also be used for many other applications, 
such as testing semi-definiteness of polynomials and copositive problems. In fact, we 
solved several difficult semi-definiteness problems efficiently by using the algorithm. Fur¬ 
thermore, applying the algorithm to copositive problems, we find an explicit expression of 
the polynomials producing open weak CADs under some conditions, which significantly 
improves the efficiency of solving copositive problems. 

The structure of this paper is as follows. In Section 2, the definitions of open weak 
delineable and open weak CAD are introduced and the problem studied in this paper 
is formally stated. In Section 3, we review basic definitions, lemmas and concepts of 
CAD. In Section 4, an algorithm for computing open weak CAD is proposed and is 
proved to be correct. In Section 5, as an application of the algorithm, we show how to 
compute open sample. In Section 6, as an application of the algorithm, a new method for 


polynomial inequality proving is proposed by combining our previous work (Han et al. 


2013) with the algorithm in Section 4. Section 7 introduces another important application 


of the algorithm to testing copositive problem. Section 8 includes several examples which 
demonstrate the effectiveness of the algorithms. 

Compared to our ISSAC’2014 paper, this paper is totally re-constructed by clarifying 
basic concepts and the problem from a new point of view. Therefore, Section 2 is totally 
re-written to start the paper by formally defining the problem studied in this paper. 
The subsequent sections are modified according to the new viewpoint. Some proofs are 
modified and added. All important new definitions and algorithms are followed by illus¬ 
trative examples. Section 7 (copositive problem) is new with substantial contents and 
new results. Section 8 contains new examples. 
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2. Problem: Open weak CAD 

In this section, we describe the problem studied in this paper. First, let us introduce 
the concept of open weak delineable. 


Definition 1 (Open weak delineable). Let S be an open set of W 1 (1 < j < n). The 
polynomial / £ R[a;i, ..., x n ] is said to be open weak delineable on S if, for any maximal 
open connected set U Cl" defined by / ^ 0, we have 

(S x R"-j) p| U ± 0 «=> Vc* £ S, (ax R n ~ j ) p U f 0. 


Let h £ R[xi,..., Xj] for j < n. We say the polynomial / is open weak delineable over h 
in R- 7 , if / is open weak delineable on any open connected component of h ^ 0 in R- 7 . 


Note that if / is analytic delineable (Collins (1975)) on S then / is also open weak 
delineable on S. However the converse is not necessarily true. 


Example 1. Let 

f = x i — ^2 + 10xl — 35^2 + 50^2 — 24 £ , aj 2 ] • 

The plot of / = 0 is given by 



Fig. 1. 

Note that / is analytically delineable and also open weak delineable on the set (—oo, —1). 
Note that / is not analytically delineable but is open weak delineable on the set (—1, oo). 
Note also that / is open weak delineable over h = x i + 1 in R. 

Definition 2 (Open weak CAD). Let / £ R[.ti, ... , x n \. A decomposition of R 7 , 1 < 
j < n, is called an open weak CAD of / in R 7 if and only if / is open weak delineable on 
every j-dirnensional open set in the decomposition. 

Example 2. Let / be the polynomial from Example 1. Let k = 1. Then 

{(- 00 ,- 1 ), [-1,-1], ( 1, oo)} 

is an open weak CAD of /. 

Finally, we are ready to state the problem precisely. 

Problem. (Open Weak CAD) Devise an algorithm with the following specification. 


3 










In: 

Out: hi, /i 2 > ■ • •, h n -1 where hj £ Z[xi,..., Xj] such that / is open weak delineable over 
hj in W. 

We call the number of open components in defined by hj ^ 0, the scale of the open 
weak CAD of / defined by hj in W. 

Remark 3. The output of the above problem is a list of “projection” polynomials, not 
an open weak CAD. However it is trivial to compute sample points in an open weak 
CAD of / from the projection polynomials. Thus, we will call the above problem “Open 
Weak CAD”. 

Example 3. Consider the following polynomial. 

In: / = ( x 3 + x 2 + x i ~ 1 )( 4 x 3 + 3x 2 + 2xi — 1) £ Z[xi,X 2 , £ 3 ] 

Out: hi = (a;i —l)(a;i + l)(29a; 2 —4a;i —24)((20a; 2 —4a;i —15) 2 + (13a; 2 —4a?i —8) 2 ) € Z[xi] 

h -2 = {%2 + x\ — l)(25x| + 12 :E 22 ;i + 20xf — 6 .T 2 — 4xi — 15) £ Z[xi,X 2 ] 

The left plot in Figure 2 below shows the open weak CAD of / produced by hi and h 2 - 
The factor (20x 2 — 4xi — 15) 2 + (13a; 2 — 4a; 1 — 8) 2 in hi does not have real root and thus 
it does not contribute to the open weak CAD. 



Fig. 2. Example 3 



Remark 4. For comparison, if we apply an open CAD algorithm on the above /, one 
would obtain the following output 

Out: hi = (xi — l)(a;i + l)(29xf — 4xi — 24){13x\ — 4xi — 8) £ Z[xi] 

/12 = (a; 2 + a; 2 — l)( 25 a:| + 12a;2a;i + 20 a; 2 — 63:2 — 4 a;i — 15 ) £ Z[a;i, X2] 
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The right plot in Figure 2 shows the open CAD of / produced by hi and /i 2 - Note that 
it has more cells than the open weak CAD (on the left). 

Remark 5. It is natural to wonder whether the multivariate discriminants of / al¬ 
ways produce open weak CADs. Unfortunately, this is not true since the discriminant 
discrim(/, [x n ,..., £j+i]) may vanish identically and thus does not always produces an 
open weak CAD of RU One may also wonder whether if the multivariate discrimi¬ 
nants of / do produce open weak CADs, then they would be the smallest open weak 
CADs. Unfortunately this is not true, either. In Example 1, it has been shown that 
X\ + 1 produces an open weak CAD of R with 2 open intervals. But the discriminant 
discrim(/, X 2 ) = —16(16a;i — 9)(xi + l) 2 produces an open weak CAD of R with 3 open 
intervals. 


3. Preliminaries 


In this section, we review some basic (known) concepts and results. If not specified, 
for a positive integer n, let x n be the list of variable (x \,..., x n ) and q„ and /3„ denote 
the point (aq,..., a n ) £ R" and (/3i,..., /3 n ) £ R", respectively. 

Definition 6. Let / £ Z[x n ], denote by lc(/, Xi) and discrim(/, aq) the leading coefficient 
and the discriminant of / with respect to (w.r.t.) Xi , respectively. The set of real zeros of 
/ is denoted by Zero(/). Denote by Zero(L) or Zero(/i,..., f m ) the common real zeros 
of L = {/ 1 ,..., fm} C Z[x n ]. The level for / is the biggest j such that deg(/, Xj ) > 0 
where deg(/, Xj) is the degree of / w.r.t. Xj. For polynomial set L C Z[x n ], I/W is the 
set of polynomials in L with level i. 


Definition 7. Let P n be the symmetric permutation group of xi,.. . , x n . Define P n q to 
be the subgroup of P n , where any element a of P n q fixes x \,... ,Xi- 1 , i.e., cr{xj ) = Xj 
for j = 1,... ,i — 1. 

Definition 8. If h£ Z[x n \ can be factorized in Z[x n ] as: 

h = ,W ■ • • ■ ■ ■ h m 2im , 

where a £ Z, t > 0 ,m > 0 , lj{i = 1, ...,i) and hi(i = 1 ,...,m) are pairwise differ¬ 
ent irreducible primitive polynomials with positive leading coefficients (under a suitable 
ordering) and positive degrees in Z[x n ], then define 

sqrfree(ft) = l\ ■ ■ ■ Iffii ■ ■ ■ h m , 
sqrfree 1 (/i) = {k,i = 1,2, ...,t}, 
sqrfree 2 (/i) = {hi, i = 1,2,, m}. 

If h is a constant, let sqrfree(Zi) = 1, sqrfree 1 (/i) = sqrfree 2 (/i) = {!}. 


In the following, we introduce some basic concepts and results of CAD. The reader 
is referred to Collins (1975); Hong ( 1990| ); Collins and Hong (1991); McCallum (1988 


1998); Brown (2001) for a detailed discussion on the properties of CAD. 
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Definition 9. (Collins, 1975 McCallum 1988) An n-variate polynomial f(x n -i,x n ) 
over the reals is said to be delineable on a subset S (usually connected) of K n_1 if (1) 
the portion of the real variety of / that lies in the cylinder SxK over S consists of the 
union of the graphs of some t > 0 continuous functions Oi < ■ ■ ■ < 8 t from S to R; and 
(2) there exist integers mi,.... m t > 1 s.t. for every a £ S, the multiplicity of the root 
8i(a) of f(a,x n ) (considered as a polynomial in x n alone) is m,. 


Definition 10. (Collins, 1975| McCallum 1988) In the above definition, the Oi are called 
the real root functions of / on S, the graphs of the Oi are called the / -sections over S, 
and the regions between successive /-sections are called f-sectors. 


Theorem 11. (McCallum, 1988, 1998) Let f(x n ,x n+ 1 ) be a polynomial in Z[® n ,a?„+i] 
of positive degree and discrim(/, x n +i) is a nonzero polynomial. Let S be a connected 
submanifold of R n on which f is degree-invariant and does not vanish identically, and 
in which discrim(/, x n +i) is order-invariant. Then f is analytic delineable on S and is 
order-invariant in each f-section over S. 

Based on this theorem, McCallum proposed the projection operator MCproj, which 
consists of the discriminant of / and all coefficients of /. 


Theorem 12. ( Brown\ 2001 1) Let f(x n ,x n +{) be a (n+l)-variate polynomial of positive 
degree in x n +i such that discrim(/, x n +i) 0. Let S be a connected submanifold ofW 1 in 
which discrim(/, x„+i) is order-invariant, the leading coefficient of f is sign-invariant, 
and such that f vanishes identically at no point in S. f is degree-invariant on S. 

Based on this theorem, Brown obtained a reduced McCallum projection in which only 
leading coefficients, discriminants and resultants appear. The Brown projection operator 
is defined as follows. 


Definition 13. (Brown 2001) Given a polynomial / G Jj[x n \, if / is with level n, the 
Brown projection operator for / is 

Bp(/, [£„]) = Res(sqrfree(/), ^^^^-,x n ). 

Otherwise Bp(/, [x n ]) = /. If L is a polynomial set with level n, then 

Bp(L, [x n ]) = (J {Res(sqrfree( /), ^^p^-,x n )} (J 


feL 


(J {Res(sqrfree(/), sqrfree(p), x n )}. 

f,g£LJ^g 


Define 


Bp(/, [x n ,x n - 1 ,... ,Xi\) 

=Bp(Bp(/, [x n , x n _i,.. - ,a: i+ i]), [aii]). 


The following definition of open CAD is essentially the GCAD introduced in Strzebonski 
(2000). For convenience, we use the terminology of open CAD in this paper. 
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Definition 14. (Open CAD) For a polynomial f(x n ) £ Z[x„\, an open CAD defined 
by f(x n ) under the order x\ -< X 2 -< • • • -< x n is a set of sample points in K” obtained 
through the following three phases: 

(1) [Projection] Use the Brown projection operator on f{x n ), let F = {/, Bp(/, [x n ]),..., 
Bp(/, [x n ,...,x 2 ])}\ 

(2) [Base] Choose one rational point in each of the open intervals defined by the real 
roots of Fb]; 

(3) [Lifting] Substitute each sample point of R* -1 for Xi_ 1 in F^ to get a univariate 
polynomial Fi(xi ) and then, by the same method as Base phase, choose sample 
points for Fi(xi). Repeat the process for i from 2 to n. 

Sometimes, we say that the polynomial set F produces the open CAD or simply, F is an 
open CAD. 


4. Algorithm: Open weak CAD 


In this section, we describe an algorithm (Algorithm 1) for computing open weak 
CAD and prove its correctness (Theorem 32). We begin by introducing a “projection” 
operator Hp. 


Definition 15 (Open weak projection operator). Let / € Z[#i, ... ,x n }. For given ?n(l < 
m < n), denote [ y] = [y 1; ..., y m \ where y.j G {a:i,..., x n j for 1 < i < m and yi 7 ^ yj for 
i 7 ^ j. For 1 < i < in, Hp(/, [y],yi) and Hp(/, [y]) are defined recursively as follows. 


B P(/. [Vi]) = Res(sqrfree(/), 


9(sqrfree(/)) 

dy% 


,Vi ) 


Hp(/, [y],Vi) = Bp(Hp (/, [y] f ), [yi]), 

Hp(/, [y]) = gcd(Hp(/, [y\, 2 / 1 ),..., Hp(/, [y],y m )), 
Hp(/, []) = /, 


where [y] i = [y u ..., yi-uVi+u ■ • •, Vm]- 


Example 4. We have 

Hp(/, [xx,xi\) = gcd (Hp (/, [xx,x 2 ], Xx ), Hp (/, [x\,x 2 ] ,x 2 )) 
Hp (/, [xi,x 2 ],x 1 ) = Bp(Hp(/, \x 2 \), [aq]) 

Hp (/, [x\,x 2 ], x 2 ) = Bp(Hp(/, [aq]), \x 2 \) 

Hp(/, [*2]) = Hp(/, [x 2 ],x 2 ) 

Hp(/, [*1]) = Hp(/, \xi],Xi) 

Hp(/, [a; 2 ],a;2) = gcd(Bp(Hp(/, [ ]), [a; 2 ])) 

Hp(/, [x\\,x\) = gcd(Bp(Hp(/, [ ]), [xi])) 

Hp(/, [ ]) = / 

Condensing the above expressions, we have 

Hp(/, [xi,x 2 ]) = gcd (Bp(Bp(/, [x 2 ]), [xi]),Bp(Bp(/, [a;i]), [aj 2 ])) 
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Algorithm 1. Open Weak CAD 
Input: f G Z[xi,... ,x n ] 

Output: hi, h 2 , ■ ■ ■, h n -i where hj G Z[xi,..., Xj] such that each hj produces an open 
weak CAD of R - 7 w.r.t. / 

1 : For all 1 < j < n and j < t < n, compute Hp(/, [x n ,..., x J+ i\,x t ) by Definition 15. 
2 : Compute hi(f), h 2 (f), ■■■, h n -i(f) where 

hj{f) = sqrfree ^ Hp(/, [x n ,..., x j+ i], x t ) 2 . 

\t=J+i / 

The Algorithm 1 (Open Weak CAD) based on the new operator Hp solves the problem 
proposed in Section 2. 

Example 5. We illustrate the Algorithm 1 using the polynomial / from Example 3. 

In: / = ( x l + x 2 + x i ~ 1 )( 4 z 3 + 3x 2 + 2#i - 1) G Z[x 1 ,x 2 ,x 3 J 

1: For all 1 < j < 3 and j < t < 3, compute Hp(/, [ 0 : 3 ,..., Xj+i], Xt) by Definition 15. 

Hp(/, [x 3 ], x 3 ) = Bp(Hp(/, [ ]), [x 3 ]) = Bp(/, [x 3 ]) 

= (x| + x\ — l)(25x| + 12 a’ 2 Xi + 20xf — 6 x 2 — 4xi — 15) 
Hp(/, [x 3 , X 2 ], X 2 ) = Bp(Hp(/, [x 3 ]), [x 2 ]) = Bp(Hp(/, [x 3 ], X 3 ), [x 2 ]) 

= (xi — l)(xi + l)(29x 2 — 4xi — 24)(13x 2 — 4xi — 8 ) 

Hp(/j [x 3 , x 2 ], x 3 ) = Bp(Hp(/, [x 2 ]), [x 3 ]) = Bp(Hp(/, [x 2 ],x 2 ), [x 3 ]) 

= (xi — l)(xi + l)(29xf — 4xi — 24)(20x 2 — 4xi — 15). 

2 : Compute h 2 (f) where 

hj(f) = sqrfree E Hp(/, [x 3 ,...,X J + i],X t ) 2 • 

\ t=J + 1 / 

h ± = sqrfree((Hp(/, [x 3 ,x 2 ], x 2 )) 2 + (Hp(/, [x 3 , x 2 ],x 3 )) 2 ) 

= (xi — l)(xi + l)(29x 2 — 4xi — 24)((20x 2 — 4xi — 15) 2 + (13xf — 4xi — 8 ) 2 ) 
h 2 = sqrfree((Hp(/, [x 3 ],x 3 )) 2 ) 

= (x| + x 2 — l)(25x 2 + 12 x 2 xi + 20xf — 6 x 2 — 4xi — 15) 

Out: hi = (xi — l)(xi + l)(29x 2 — 4xi — 24)((20x 2 — 4xi — 15) 2 + (13xf — 4xi — 8 ) 2 ) 
h 2 = {x\ + x\ — l)(25x| + 12 x 2 xi + 20xf — 6 x 2 — 4xi — 15) 

Remark 16. Although Hj is more complicated than any Hp(/, [x n ,..., Xj+i], Xt), 

Zero (Hj) C Zero(Hp(/, [x n ,..., Xj+i], x j+ i)) C Zero(Bp(/, [x n ,... ,x j+ i])). 

That means, for every open cell C of open CAD produced by Brown’s operator, there 
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exists an open cell C of open weak CAD produced by Hj such that C' C C . Thus, the 
scale of open weak CAD is no bigger than that of open CAD. 

Remark 17. In Algorithm 1, the scale of the open weak CAD of / defined by hj in 
R 7 is not always the smallest. For example, let / be the polynomial in Example 1, then 
hi = (16a;i — 9)(si + 1 ), and / is open weak delineable over x\ + 1 , as mentioned earlier. 

We now prove the correctness of Algorithm 1. 

Definition 18. (Open sample) A set of sample points Sf C R fe \ Zero(/) is said to be 
an open sample defined by f{xk) £ Z[x^] in R fc if it has the following property: for every 
open connected set U C R fc defined by / ^ 0, Sf fl U ^ 0. 

Suppose g(xk) is another polynomial. If Sf is an open sample defined by f(xk) in R fc 
such that g(a) ^ 0 for any a £ Sf, then we denote the open sample by S/, g . 

As a corollary of Theorems 11 and 12, a property of open CAD (or GCAD) is that 
at least one sample point can be taken from every highest dimensional cell via the open 
CAD (or GCAD) lifting phase. So, an open CAD is indeed an open sample. 

Obviously, there are various efficient ways to compute S/ jS for two given univariate 
polynomials f,g £ Z[x\. For example, we may choose one rational point from every open 
interval defined by the real zeros of / such that g does not vanish at this point. Therefore, 
we only describe the specification of such algorithms SPOne here and omit the details of 
the algorithms. 

Definition 19. Let a.j = (au,..., atj) £ R - 7 and S' C K be a finite set, define 

atj EE S = {(ai, ...,aj,ft) \ ft £ S}. 


Algorithm 2. SPOne 

Input: Two univariate polynomials f,g £ Z[x\ 

Output: Sf t g, an open sample defined by f(x) in R such that g{a) ^ 0 for any a £ Sf, g 


Algorithm 3. OpenSP 

Input: Two lists of polynomials Li = [f n {x n ),.. L 2 = [ g n {x n ), • • .,gj(xj)], and 

a set of points S in W 
Output: A set of sample points in R” 

1: 0:=S 

2 : for i from j + 1 to n do 

3: P := 0 

4: for a in O do 

5 : p := P{J{a ffl SP0ne(/j(a, Xj),gi(a, £,))) 

6: end for 

7: 0:=P 

8 : end for 
9: return O 
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Remark 20. For a polynomial f(x n ) G Z[x n ], let 

Bi = [/,Bp(/, [x„]),.... Bp(/, [x n ,..., x 2 ])], 

S = SP0ne(Bp(/, [x n ,... ,x 2 ]), 1), 

then OpenSP(i?i, B 2 , S) is an open CAD (an open sample) defined by f(x n ). 

We will provide in this section a method which computes two lists C 1 and C 2 where 
the polynomials in C\ are factors of corresponding polynomials in B 1 and will prove that 
any OpenSP(Ci, C 2 , Sf ) is an open sample ofK” defined by f(x n ) for any open sample 
Sf^g o in R 7 where fj G C\ and gj G C 2 . 

Remark 21. The output of OpenSP(Li, L 21 S) is dependent on the method of choosing 
sample points in Algorithm SPOne. In the following, when we use the terminology “any 
QpenSP(L!, L- 2 , S )”, we mean “no matter which method is used in Algorithm SPOne for 
choosing sample points”. 

Definition 22. (Open delineable) Let 

Ll [f n {x n 'j, fn—l(x n —l), ■ ■ ■ , fj(Xj)\, (1) 

L 2 = \9n{Xn)i Qn— 1 (x n — i),. . . , gj {Xj )] ( 2 ) 

be two polynomial lists, S an open set of R s (s < j) and S' = S x R J_S . The polynomial 
f n (x n ) is said to be open delineable on S w.r.t. L\ and L 2 , if ^ 0 for any maximal 

open connected set U C M” defined by /„ 7 ^ 0 with [/f'|(S' / x K” - - 5 ) ^ 0 and any 
A = OpenSP(Li, L 2 , {<*}) where a. G S' is any point such that fj(at)gj(a) ^ 0. 

Remark 23. Let s = j in Definition 22, it could be shown that if f n {x n ) is open 

delineable on S' w.r.t. L\ and L 21 then f n (x n ) is open weak delineable on S'\Zero{fjgj}. 

Suppose f n (x n ) is a squarefree polynomial in Z[x n ] of positive degree and S C K n_1 
is an open connected set in which Bp(/„, \x n \) is sign-invariant. According to Theorem 
11 and Theorem 12, f n is analytic delineable on S. It is easy to see that f n is open 
delineable on S w.r.t. [/ n ,Bp(/„, [*„])] and [/„,Bp(/„, [*„])]. 

Open delineability has the following four properties. 

Proposition 24. (open sample property) Let Li,L 2 be as in Definition 22. If f n (x n ) is 
open delineable on every open connected set of fj{xj) 0 w.r.t. L\ and L 2 , then for any 
open sample Sf. 9j in K- 7 , any A = 0 penSP(Li, L 2 , Sf. g A is an open sample defined by 
f„(x n ) in K". 

Proof. For any open connected set U C K” defined by f n ^ 0, there exists at least 
one open connected set S C defined by fj 7 ^ 0 such that [/(^(S 1 x R” - - 7 ) 7 ^ 0. 
Since /„ is open delineable on S w.r.t. L\ and L 2 , we have Af)U 7 ^ 0 for any A = 
OpenSPf L 1 ,L 2 ,S f ^g :j ). □ 

Proposition 25. (transitive property ) Let L\, L 2 , S, S' be as in Definition 22. Suppose 
that there exists k(j < k < n) such that fk(xk) is open delineable on S w.r.t. L'{ = 
[fk(xk),---,fj(xj)] and L'f = [g k (x k ),..., g j (x j )}, and f„(x n ) is open delineable on 
every open connected set of f k (x k ) 7 ^ 0 w.r.t. L\ = [f n (x n ),..., f k {x k )] and L' 2 = 
[g n (x n ), ■ ■ ■, gk(x k )\. Then f n (x n ) is open delineable on S w.r.t. L\ and L 2 . 
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Proof. Let a £ S' be any point such that fj(a)gj(a ) 7 ^ 0, for any A = OpenSPfXi, L 2 , {a}), 
we have A = 0penSP(L , 1 , L 2 ,A') where A! = 0penSP(L", L 2 , {«})■ For any open con¬ 
nected set U C K™ defined by f n 7 ^ 0 with U p|(S" x K " - - 7 ) 7 ^ 0, there exists an open 
connected set V C defined by fk 0 with U f)(V xR n ~ k ) 7 ^ 0 and V f'|(<S'xii fc-s ) 7 ^ 0. 
Now we have A' 7 ^ 0 since fk(xk) is open delineable on S w.r.t. L'{ and L 2 . And 
then, Af] U 0 is implied by x K ra_fc ) 7 ^ 0 since f n {x n ) is open delineable on V 

w.r.t. L'| and Lj. n 

Proposition 26. (nonempty intersection property) Let L\, L 2 be as in Definition 22. 
For two open sets S 1 and S 2 o/K s (s < j) with Sif)S 2 ^ 0, if f n (x n ) is open delineable 
on both S 1 and S 2 w.r.t. Li and L 2 , then f n (x n ) is open delineable on Si (J S 2 w.r.t. L\ 
and L 2 • 

Proof. For any a 1 € £ 1 , a 2 € £ 2 , a 3 € Si with fj(oti)gj(oti) 7 ^ 0, any At = 
QpenSP(Li, L 2 , {o:*})) and open connected set U C K ra defined by f n 7 ^ 0 , we have 
Uf){Si xR n ~ s ) ^0) ^ Aif)U AzC\U A 2 f)U <=>U0(S 2 x 

R”- s ) 7^ 0. □ 

Proposition 27. (union property) Let Li,L 2 be as in Definition 22. For a £ P n ,j+ 1 , 
denote y n = (yi, ■.. ,y n ) = cr(x n ) and y l = (yi ,..., yf). Let L[ = [f n (x n ), p n -i(y n -i), 
..., Pj(yj)] and L' 2 = [q n (x n ), q n -i(y n -l), ■■■■■ QjiVj)] where Pifa) and q^yf) are 
polynomials in i variables. 

For two open sets S 1 and S 2 ofW, if (a) f n {x n ) is open delineable on both S\ and 
S 2 w.r.t. L\ and L 2 , (b) f n (x n ) is open delineable on S 1 [JS 2 w.r.t. L[ and L' 2 , and (c) 
Pj(Vj)Qj(Vj) vanishes at no points in Si [J S 2 , then f n {x n ) is open delineable on Si (J S 2 
w.r.t. L 1 and L 2 . 

Proof. Let a 1 £ S 1 , 01.2 €E S 2 be two points such that gjPjqj(a t ) 7 ^ 0 for t = 1,2. Let 
At = OpenSPfXi, L 2 , {ct t }) and A't = 0 penSP(L , 1 , L 2 , {a t }). For any open connected set 
U defined by /„ 7 ^ 0 with U [){oli x R” - - 7 ) 7 ^ 0, then Aif)U 7 ^ 0 and A' if)U 7 ^ 0. 
Since f n {x n ) is open delineable on S 1 (J S 2 w.r.t. L\ and L' 2l we have A' 2 p| U 7 ^ 0 which 
implies that U P|(5 2 x JR” - - 7 ) 7 ^ 0 and A 2 H U 7 ^ 0. Therefore, f n (x n ) is open delineable 
on (Si (J S' 2 )\Zero(pj< 7 . ) ') w.r.t. Li and L 2 . Since Pjqj does not vanish at any point of 
Si (J S 2 , fn{x n ) is open delineable on Si (J S 2 w.r.t. L\ and L 2 . □ 

Example 6. We illustrate Proposition 27 using the polynomial / from Example 3. 

f = (xl+xl+xf- 1)(4x 3 + 3x 2 + 2xi — 1) e Z[xi,x 2 , £ 3 ] 

Let 


F'l [f 1 fx 3 , fx 2 X 2] 
L'l = [f,fx 2 Jx 2 x 3] 

L 2 = Li 

T' — T' 

^ 2 — ^1 


11 



where 


f X3 = (x% + x\ — l)(25x| + 12 x 2 Xi + 20x} — 6 x 2 — 4xi — 15), 
fx 3x2 = ( x \ - 1)0*4 + l)(29xf - 4x x - 24)(13x? - 4xi - 8 ), 
f X2 = (x§ + x\ — l)(25x§ + 16 x 3 X 1 + 13xi — 8 x 3 — 4xi — 8 ), 
fx 2 x 3 = (®i - l)(xi + l)(29x? - 4xi - 24)(20x? - 4xi - 15). 

Let 

2 6V3 2 6^3^ 

13 ~ ITT’ 13 + ~W) 

2 6%/3 1 

13 + "lT’ To + ~5~ 

be two open intervals in xi-axis, where Xi = ^ is one of the real roots of the 

equation f x3X2 = 0. By typical CAD methods, Si and S 2 are two different cells in xi- 
axis. 

It could be deduced easily by Theorem 11, Theorem 12, and Proposition 25 that the 
conditions (a) and (b) of Proposition 27 are satisfied. Since f X 2 x 3 vanishes at no points 
in S 1 \JS 2 , condition (c) is also satisfied. 

By Proposition 27, / is open delineable on Si IJ S 2 w.r.t. Li and L 2 . Roughly speaking, 
the real root of Xi — ^ would not affect the open delineability, thus we could combine 
the two cells Si and S 2 . 

For the same reason, the real roots of 

(13xi - 4xi - 8)(20xi - 4xi - 15) 
would not affect the open delineability either. 




Lemma 28. (Han et al. 
open set U of] 


2013) Let f and g be coprime in Z[x n J. For any connected 


the open set V = U\Zero(f 1 g) is also connected. 


Lemma 29. Let f = gcd(/i,..., f m ) where fi £ Z[x n ], i = 1,2,, m. Suppose f has 
no real zeros in a connected open set U C then the open set V = /7\Zero(/i,..., f m ) 
is also connected. 


Proof. Without loss of generality, we can assume that / = 1. If m = 1, the result 
is obvious. The result of case m = 2 is just the claim of Lemma 28. For m > 3, let 
9 = gcd(/i,..., f m -i) and = f t /g (i = 1 1), then gcd (f m ,g) = 1 and 

gcd(gi,.. .,g m - 1 ) = 1- Let A = Zero(/i,... ,/ m ), B = Zero(gi,. .., g m -i) U Zero^, / m ). 
Since A C B, we have U\B C U\A. Notice that the closure of U\B equals the closure 
of U\A, it suffices to prove that U\B is connected, which follows directly from Lemma 
28 and induction. □ 


Definition 30. Define 

Hp(/, i) = {/, Hp(/, [x„]),..., Hp(/, [x n ,..., x^)}, 

Hp(/, i) = {/, Hp(/, [x„], X n ), . . . , Hp(/, [x„, ... , Xi], Xi)}. 
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As a Corollary of Theorem 11 and Theorem 12, we have 


Proposition 31. Let / € Z[x n \ be a squarefree polynomial with level n. Then f is open 
delineable on every open connected set defined by Bp(/, [x n ]) 0 in R " _1 w.r.t. Hp(/, n) 

and Hp (/, n). 

The following Theorem is the main result of this paper, which guarantees the correct¬ 
ness of Algorithm 1. 

Theorem 32. Let j be an integer and 2 < j < n. For any given polynomial f(x n ) € 
Z [x n \ and any open connected set U C o/Hp(/, [ x n ,..., Xj]) 7 ^ 0, let S = t/\Zero( 
{Hp(/, [ x n ,... ,Xj],Xt) | t = j,... ,?r}). Then f{x n ) is open delineable on the open con¬ 
nected set S w.r.t. Hp(/, j) and Hp(/, j). As a residt, f is open weak delineable over 
in R-? , and Algorithm 1 is correct. 

Proof. First, by Lemma 29, S is open connected. We prove the theorem by induction on 
k = n — j. When k = 0, it is obviously true from Proposition 31. Suppose the theorem is 
true for all polynomials g(xk) € Z[x;,.] with k = 0,1,..., n—i—1. We now consider the case 
k = n — i. Let [z] = [x„,..., xf\. For any given polynomial f(x n ) G Z[cc„], let U C R I_1 
be an open connected set of Hp(/, [z]) ^ 0 and S = U\Zero({Hp(/, [z],x t ) | < = n}). 

For any point at £ S with Hp(/, [z\,Xi)(at) ^ 0, there exists an open connected set 
S a C R ' i_1 such that at G S a and 0 ^ Hp(/, [z], Xi)(S a ). By induction, Hp(/, [ x n ,..., x <+iD 
is open delineable on S a w.r.t. {Hp(/, [z])} and {Hp(/, [z],£j)}. By induction again and 
the transitive property of open delineable (Proposition 25), / is open delineable on S a 
w.r.t. Hp(/,«) and Hp(/,i). 

For any point at G S with Hp(/, [z], x*)(a) =0, there exists an i' such that n > i' > i+1 
and Hp(/, [z\,Xi’){at) ^ 0. Thus there exists an open connected set S' a of R I_1 such that 
a G S' a and 0 ^ Hp(/, [z], xp)(S' a ). Let a G P n ^ with a{xi) = Xi <, in such case, f(a(x n )) 
is open delineable on S' a w.r.t. Hp(/(cr(a: n )), *) and Hp(/(cr(a:„)), *). For any /3 G S' a with 
Hp(/, [z], Xi)(/3) ^ 0, there exists an open connected set S'" C S' a and / is open delineable 
on S" w.r.t. Hp (f,i) and Hp(/, i). From union property of open delineable (Proposition 
27), / is open delineable on S' a w.r.t. Hp(/, i) and Hp(/,«). 

To summarize, the above discussion shows that for any point at G S, there exists 
an open connected set S a C S such that at G S a and / is open delineable on S a 
w.r.t. Hp(/, i) and Hp(/, *). By the nonempty intersection property of open delineable 
(Proposition 26) and the fact that S is connected, f(x n ) is open delineable on S w.r.t. 
Hp(/, i) and Hp(/, i) as desired. 

Therefore, the theorem is proved by induction. The last statement of the theorem 
follows from the fact Zero({Hp(/, [ x n ,..., Xj\,x t ) \ t = j ,..., n}) = Zero(77j_i(/)). □ 


5. Application: Open Sample 

As a first application of Theorem 32, we show how to compute open sample based on 
Algorithm 1. 

Definition 33. A reduced open CAD of f{x n ) w.r.t. [x n ,... ,Xj+ 1 ] is a set of sample 
points in R" obtained through Algorithm 4. 
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Algorithm 4. ReducedOpenCAD 

Input: f(x. n ) G Z[x n \ and an open sample [*»,->^ + 1 ]) >H p(/,[*n,.-,*^'1.^+0 in 

Output: A set of sample points in K” 

1: O := ^Hp(/,[a; n ,...,Xj+i]),Hp(/,[x n ,...,Xj.fi] , 0 ^+ 1 ) 

2 : for i from j + 2 to n + 1 do 

3: P := 0 

4: for a in O do 

5: if i <n then 

6: P := PU(aB3S p 0ne(Hp(/, [x n , ■ ■ ., *»])(<*, Xi-i), Hp(/, [x„,.. Xi)(a, ®i_i))) 

7: else 

8: P := PU( a EB SP0ne(/(a, x n ), f(a,x„))) 

9: end if 

10 : end for 

11 : O - P 

12 : end for 
13: return O 

The following corollary of Theorem 32 shows that the reduced open CAD owns the 
property of open delineability. 

Corollary 34. A reduced open CAD of f(x n ) w.r.t. [x n , ... ,Xj+ 1 ] is an open sample 
defined by f(x n ). 

Example 7. We illustrate the main steps of Algorithm 4. using the polynomial / from 
Examples 3 and 5. 

In: / = (x§ + x\ + x\ — 1)(4x3 + 3x 2 + 2#i — 1) G Z[xi,x 2 ,x 3 ] 

^ , Hp(/,[x 3 ,X 2 ]),Hp(/,[x 3 ,X 2 ],X 2 ) ~ 2, —32,0, 54 , 2 } in M. 

1: O := {—2, —1|, 0, ||, 2} (O has 5 elements, ap,..., 05 ) 

3: P := 0 

6: P := PU(ar ffl SP0ne(Hp(/, [x 3 ])(c*i,x 2 ),Hp(/, [x 3 ],x 3 )(ai,x 2 ))) 

P :=P|J(« 2 ffl SP0ne(Hp(/, [x 3 ])(a 2 ,x 2 ),Hp(/, [x 3 ],x 3 )(a 2 , x 2 ))) 

P ■= ^U ( q 5 ffl SP0ne(Hp(/, [x 3 ])(a 5 , x 2 ),Hp(/, [x 3 ],x 3 )(a 5 ,x 2 ))) 

11 O := P (O now has 13 elements, ai 3 ) 

3: P := 0 

8: P ■= P U(ar ffl SP0ne(Hp(/)(a 1 , x 3 ), Hp(/)(o: 1 , x 3 ))) 

P := PU («2 ffl SPDne(Hp(/)(a 2 ,x 3 ),Hp(/)(a 2) x 3 ))) 

P := PU( a i3 ffl SP0ne(Hp(/)(ai3,x 3 ),Hp(/)(ai 3 ,x 3 ))) 

11: O := P (O now has 36 elements, cti,..., a 36 ) 
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Out O 


Remark 35. As an application of Theorem 32, we could design a CAD-like method to 
get an open sample defined by f{x n ) for a given polynomial f(x n ). Roughly speaking, 
if we have already got an open sample defined by Hp(/, [x n ,... ,Xj]) in R J , according 
to Theorem 32, we could obtain an open sample defined by / in M". That process could 
be done recursively. 

In the definition of Hp, we first choose m variables from {xi, ..., x n }, compute all 
projection polynomials under all possible orders of those m variables, and then compute 
the gcd of all those projection polynomials. Therefore, Theorem 32 provides us many ways 
for designing various algorithms for computing open samples. For example, we may set 
m = 2 and choose [x n ,x n -i], [x n -2,x n -3,], etc. successively in each step. Because there 
are only two different orders for two variables, we compute the gcd of two projection 
polynomials under the two orders in each step. Algorithm 5 is based on this choice. 


Algorithm 5. HpTwo 

Input: A polynomial / £ Z[x n \ of level n. 

Output: An open sample defined by /, i.e., a set of sample points which contains at least 
one point from each connected component of / ^ 0 in R" 


1 

9 

— /; Li {}; L 2 

:= 

{}; 

2 

while i > 3 do 



3 


L\ := Li |J Hp(<7: 

i — 

i); 

4 


L 2 := L 2 \Jttp{g, 

i — 

i); 

5 


9 ■= [xi,Xi- 

-i]) 


6 


i := i — 2; 



7 

end while 



8 

if 

i = 2 then 



9 


Li := Li |J Hp(<7, 

*); 


10 


L 2 ■= T 2 U H P(5: 

*); 


11 


9 ■= Hp(s, N); 



12 

end if 



13 

S: 

=SP0ne(4 1 J ,4 1 ) 



14 

C 

= OpenSP(Li, L 2 , 

sy, 


15 

return C. 




Remark 36. If Hp(/, [x n , x n -\)] ^ Bp(/, [x„, x n -i\) and n > 3, it is obvious that the 
scale of projection in Algorithm 5 is smaller than that of open CAD in Definition 14. 


Remark 37. It should be mentioned that there are some non-CAD methods for com¬ 
puting sample points in senri-algebraic sets, such as critical point method. For related 


work, see for example, Basu et al. (19981; Safey El Din and Schost (2003); Hong and 
Safey El Din (2012). Mohab Safey el Din developed a Maple package RAGlitp | based o~n 


their work, which can test semi-definiteness of polynomials. 


1 http://www-polsys.lip6.fr/ safey/RAGLib/distrib.html 
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6. Application: Polynomial Inequality Proving 


In this section, we combined the idea of Hp and the simplified CAD projection operator 
Np we introduced previously in Han et al. (2013), to get a new algorithm for testing semi¬ 
definiteness of polynomials. 


Definition 38. (Han et ah, 2013) Suppose / £ Z[x n ] is a polynomial of level n. Define 


Oc(/, x n ) = sqrfree 1 (lc(/, x n )), Od(/, x„) = sqrfree 1 (discrim(/, x n )), 
Ec(/, x n ) = sqrfree 2 (lc(/,x n )),Ed(/,x n ) = sqrfree 2 (discrim(/, x„)), 
Ocd(/, X n ) = Oc (/, X n ) U Od(/, X n ), 

Ecd(/, x n ) = Ec(/, Xn) U Ed(/, Xn). 


The secondary and principal parts of the projection operator Np are defined as 

NPi(/, K]) =Ocd(/, x n ), 

%(/> w) ={ n g }■ 

geEcd(/,a:„)\Ocd(/,a;„) 

If L is a set of polynomials of level n, define 

Npi(L, [x„]) = (J Ocd(y, x n ), 

geL 

Np 2 (i, W) = IJ { II M- 

g£L /iGEcd(^,a: Tl )\Np 1 (L,[x Tl ]) 


Based on the projection operator Np, we proposed an algorithm, DPS, in (Han et al. 


2013) for testing semi-definiteness of polynomials. Algorithm DPS takes a polynomial 
f(x n ) £ Z\x n \ as input, and returns whether or not /( x n ) > 0 on K". The readers are 
referred to (I Ian et al. 2013) for the details of DPS. 

The projection operator Np is extended and defined in the next definition. 


Definition 39. Let / £ Z[x\,... ,x n \ with level n. Denote [y\ = [y±,... ,y m ], for 1 < 
m < n, where yi £ {aq,..., x n } for 1 < i < m and yi ^ yj for i ^ j. Define 

Np(/, [Xi]) = Np 2 (/, [a.’;]), Np(/, [xi], Xi) = 9- 

SreNpi(/,[rci]) 

For m(m > 2) and i( 1 < i < m), Np(/, [y],yi) and Np(/, [y]) are defined recursively as 
follows. 

Np(/, [ y],Vi ) = Bp(Np(/, [y) t ),yi), 

Np(/, [y]) = gcd(Np(/, [y], yi),..., Np(/, [y\, y m )), 

where [y\ i = [y 1; ..., yi-i, y i+1 , ...,y m ). Define 

Np(/i i) = {/, Np(/, [x n ]), . . . , Np(/, [x n ,. .. , Xi)]}, 

and 

Np(/, i) = {/, Np(/, [x ra ], X n ), . . . , Np(/, [x n , • • • , Xi],Xi)}. 
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Theorem 40. (Han et al. 2013) Given a positive integer n > 2. Let f £ Z[x n ] be 
a non-zero squarefree polynomial and U a connected component of Np (/, [s ra ]) 7 ^ 0 in 
M n_1 . If the polynomials in Np 1 (/, [a; n ]) are semi-definite on U, then f is delineable on 

V = U\\Jh& ?1 (/,[*„]) Zero (h). 


Lemma 41. (Han et al., 2013) Given a positive integer n > 2. Let f £ Z[x n ] be 


squarefree polynomial with level n and U a connected open set of Np(/, [x n \) 7 ^ 0 in 
R n_1 . If f(x. r j) is semi-definite on U x R, then the polynomials in Np 1 (/, [x n \) are all 
semi-definite on U. 


Now, we can rewrite Theorem 40 in another way. 

Proposition 42. Let f £ Z[x n ] be a squarefree polynomial with level n and U a connected 
component o/Np(/, [x n ]) 7 ^ 0 in If the polynomials in Np 1 (/, [&„]) are semi-definite 

on U, then f is open delineable on U w.r.t. Np(/, n) and Np (f,n). 

Theorem 43. Let j be an integer and 2 < j < n. For any given polynomial f(x n ) £ 
Z[x n ], and any open connected set U o/Np(/, [x n ,..., Xj\) 7 ^ 0 in R J , let 

S = £/\Zero({Np(/, [x n ,...,Xj\,x t ) \ t=j,...,n}). 

If the polynomials in NpiCy, [ 0 ;^—^]) are all semi-definite on U x f(x n ) is 

open delineable on S w.r.t. Np (f,j) andHp(f,j). 


Proof. Notice that the proof of Theorem 32 only uses the properties of open delineable 
(Propositions 24-27) and Proposition 31(/ is open delineable on every open connected set 
defined by Bp(/, [x„]) 7 ^ 0 in M " -1 w.r.t. Hp(/, n) and Hp(/, n)). According to Proposition 
42, / is open delineable on every open connected set defined by Np(/, [x„]) 7 ^ 0 in K n_1 
w.r.t. Np(/, n) and Np(/, n). The same proof of Theorem 32 will yield the conclusion of 
the theorem. □ 


Theorem 43 and Proposition 41 provide us a new way to decide the non-negativity of 
a polynomial as stated in the next theorem. 

Theorem 44. Let f £ Z[x n ] be a squarefree polynomial with level n and U a connected 
open set o/Np(/, [. x n ,..., Xj]) 7 ^ 0 in R- 5 . Denote S = C/\Zero({Np(/, [x n , ... ,Xj],Xt) \ 

t = j,...,n}). The necessary and sufficient condition for f(x n ) to be positive semi- 
definite on U x K™ _ 1 +1 is the following two conditions hold. 

(1) The polynomials in [J( = q Np 1 (/, [x n -i\) are all semi-definite on U x R n_ h 

(2 ) There exists a point a £ S such that f(a,Xj,...,x n ) is positive semi-definite on 
R n ~i +1 . 


Based on the above theorems, it is easy to design some different algorithms (depending 
on the choice of j) to prove polynomial inequality. For example, the algorithm PSD-HpTwo 
for deciding whether a polynomial is positive semi-definite, which we will introduce later, 
is based on Theorem 44 when j = n — 1 (Proposition 45). 
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Proposition 45. Given a positive integer n > 3. Let / £ Z[x n \ be a squarefree polyno¬ 
mial with level n and U a connected open set of Np(/, [x n ,x n -i\) ^ 0 in R n ~ 2 . Denote 
S = Z7\Zero(Np(/, [x n ,x n -i\,x n ),Np{f, [x n ,a:„_i],a: n _i)). 

The necessary and sufficient condition for /( x n ) to be positive semi-definite on U x R 2 
is the following two conditions hold. 

(1) The polynomials in either Np 1 (/, [i„]) or Np 1 (/, [x n -i]) are semi-definite on U xl. 

(2 ) There exists a point a £ S such that f(a 1 x n -i,x n ) is positive semi-definite on R 2 . 


Algorithm 6. PSD-HpTwo 

Input: An irreducible polynomial / £ T\x n \. 

Output: Whether or not Vcc n £ R”, f(a . n ) > 0. 

1: if n < 2 then 

2: if DPS(/(a; n ))=false then return false 

3: end if 

4: else 

5: Li :=N Pl (/, [a„])U N Pl(/. fcn-l]) 

6 : L 2 := Np(/, [x n ,x n -i]) 

7: for g in L\ do 

8: if PSD-HpTwo(g) =false then return false 

9: end if 

10: end for 

11 : C n - 2 := A reduced open CAD of L 2 w.r.t. [x n - 2: ... ,x 2 ], which satisfies that 

Zero(Np(/, [x n ,x n -i\,x n ),Np(f, [x n ,x n -i),x n -i)) nC '„_ 2 =0. 

12 : if 3 q :, 1 _2 £ C n - 2 such that DPS(/(q: T i_ 2 ,x n _i, a; n ))=false then return false 

13: end if 

14: end if 

15: return true 


7. Application: Copositive problem 

Definition 46. A real n x n matrix A n is said to be copositive if > 0 for every 

nonnegative vector x n . For convenience, we also say the form x n A n x^ is copositive if 
A n is copositive. 


The collection of all copositive matrices is a proper cone; it includes as a subset the 
collection of real positive-definite matrices. For example, xy is copositive but it is not 
positive semi-definite. 

In general, to check whether a given integer square matrix is not copositive, is NP- 


complete (Murty and Kabadi 1987). This means that every algorithm that solves the 


problem, in the worst case, will require at least an exponential number of operations, 
unless P=NP. For that reason, it is still valuable for the existence of so many incomplete 
algorithms discussing some special kinds of matrices (Parrilo 2000). For small values 


of n (< 6 ), some necessary and sufficient conditions have been constructed (Andersson 


et al.|[l995| Hadeler 1983). We refer the reader to (Hiriart-Urruty and Seeger 20101 for 


a more detailed introduction to copositive matrices. 
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From another viewpoint, this is a typical real quantifier elimination problem, which can 
be solved by standard tools of real quantifier elimination ( e.g ., using typical CAD). Thus, 
any CAD based QE algorithm can serve as a complete algorithm for deciding copositive 
matrices theoretically. Unfortunately, such algorithm is not efficient in practice for the 
CAD algorithm is of doubly exponential time complexity. 

To test the copositivity of the form x n A n x is equivalent to test the nonnegativity 
of the form (x 2 , ..., x 2 )A n (xf ,..., x 2 ) T . In this section, we give a singly exponential in¬ 
complete algorithm with time complexity 0 (n 2 4") based on the new projection operator 
proposed in the last section and Theorem 44. 

Let us take an example to illustrate our idea. Let 

F := ax 4 + bx 2 y 2 + cy 4 + dx 2 + ey 2 + /, 

be a squarefree polynomial, where a, b, c € Z, d,e, f £ Z [z n \ and a 7 ^ 0, c ^ 0. 

To test the nonnegativity of F, we could apply typical CAD-based methods directly, 
i.e., we can use Brown’s projection operator. In general, we have 

Bp(F, [x]) = (cy 4 + ey 2 + f)a(4acy 4 + 4aey 2 + 4af — b 2 y 4 — 2 by 2 d — d 2 ). 

We then eliminate y, 

Bp (F, [x,y]) = Bp(Bp(F, [x]),y) 

=fac(4fc — e 2 )(d 2 c — edb + fb 2 )(4af — d 2 )(4ac — b 2 )(4afc — ae 2 — d 2 c + edb — fb 2 ). 

If d, e are polynomials of degree 2 and / is a polynomial of degree 4 (copositive problem 
is in this case), the degree of the polynomial Bp(F, [x, y}) is 20 while the original problem 
is of degree 4 only. That could help us understand why typical CAD-based methods do 
not work for copositive problems with more than 5 variables in practice. 

Now, we apply our new projection operator. Notice that 

3F 

Res (F, — ,x) = 16 (cy 4 + ey 2 + f)F 2 , 
ox 

where F\ = a(4acy 4 + 4 aey 2 + 4 af — b 2 y 4 — 2 by 2 d — d 2 ). 

If cy 4 + ey 2 + f and F\ are nonzero and squarefree, Np(F, [x]) = F\. Thus, in order to 
test the nonnegativity of F, we only need to test the semi-definitness of cy 4 + ey 2 + /, 
choose sample points defined by Np(F, [x]) 7 ^ 0 (we also require that cy 4 +ey 2 + f does not 
vanish at those sample points) and test the nonnegativity of F at these sample points. 
On the other side, 

Res (F, = 16(ax 4 + dx 2 + f)F 2 , 

where F 2 = c(4cx 4 a + 4cdx 2 + 4/c — b 2 x 4 — 2bx 2 e — e 2 ). 

Similarly, if ax 4 + dx 2 + f and F 2 are nonzero and squarefree, Np(F, [y]) = F 2 . In order 
to test the nonnegativity of F, we only need to test the semi-definitness of ax 4 + dx 2 + /, 
choose sample points defined by Np(F, [y]) 7 ^ 0 (we also require that ax 4 +dx 2 + f does not 
vanish at those sample points) and test the nonnegativity of F at these sample points. 

Under some “generic” conditions (i.e., some polynomials are nonzero and squarefree), 
we only need to test the semi-definitness of ax 4 + dx 2 + / and cy 4 + ey 2 + /, choose 
sample points T 2 defined by Np(F, [x,y]) = gcd(Bp(Np(F, [x]), y), Bp(Np(F, [y]), x)) = 
(4ac — b 2 )(4afc — ae 2 — d 2 c + edb — fb 2 ) 7^ 0 (we also require that Res(Np(F, [x]), y) does 
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not vanish at T 2 ), obtain sample points T\ defined by Np (F, [a;]) ^ 0 from T 2 (we also 
require that cy 4 + ey 2 + / does not vanish at Xj) and test the nonnegativity of F at Tj. 

Again, if d, e are polynomials of degree 2 and f is a polynomial of degree 4, both the 
degree of Np(F, [a;]) and Np(F, \x,y\) are exactly 4. It indicates that our new projection 
operator may control the degrees of polynomials in projection sets. Moreover, we point 
out that 

Np(F, [a:, y]) = (4 afc — ae 2 — d 2 c + edb — fb 2 ) 

' ( 

= 4det 

V 

Before giving the result, we introduce some new notations and lemmas for convenience. 


b d 

\ 

a 2 2 


6 „ e 


2 c 2 


d e £ 

2 2 J 

/ 


Definition 47 (Sub-sequence). An array I is called a sub-sequence of sequence {1, ..., n} 
if for any i-th component of I, I[i] € {!,..., n} and I[i] < I[i + 1] for i = 1,..., |/| — 1. 


For a sub-sequence / of {1,..., n} with m = |/|, denote xj = [aqm,..., a xj = 
{xi | i ^ 1} and Aj = (a/m ,/m )i<i,j<m the sub-matrix of A n . 

Let 

n 

y*(*^n) — ^ [ (M., 3 X%X :i “t" ^ ^ (ai,n+1 "t ,i)Xi -f l,n-\-l — (Xn , 1 )A n -|_i (a? n , 1) (3) 

be a quadratic polynomial in x n , where A n+ 1 = (a,i,j)t n +i)x(n-\-i), ciij = a^i € Z[z s ] for 
1 < i,j < n + 1. Set F(x n ) = f(xf, ..., a;^). It is not hard to see that (please refer to the 
proof of Theorem 49), for a given sub-sequence I of {1,..., n} with length in, there exist 
some polynomials p±,... ,p m + 1 G Z[z s ,xj] such that F{x n ) = ( x 2 , 1 )Pi(x 2 , 1) T where 


Pi 


A/ ( Pi,---,p m ) T 
(Pl> ‘ ‘ ‘ i Pm) Pm+1 


For convenience, we denote Pn_ )m i by P m +i(a: m +i,... ,x n ) or simply P m+ i. In partic¬ 
ular, P „+1 = A n+ 1 and F(x n ) = (xj,.. .,x 2 , l)P. i+1 (x\,.. .,x 2 , 1) T . 


Example 8. Suppose F(aq, x 2 , 3 : 3 ) = xf+2x2+4xfx% — 2xfx§+4x%x§ + 8xfz 2 + 5x3+z 4 . 
If I = [ 1 , 2 ], then */ = [xi,x 2 ],W[ = {a; 3 }, 


Pi 


1 2 4 z 2 — x 2 

2 2 2x\ 

4z 2 — x\ 2x 2 5a;| + 2 4 


If I = [1], then x / = [aq],cc/ = {a; 2 ,a: 3 }, 


P/ 


1 2x\ + 4 z 2 — x 2 

2x 2 + 4z 2 — x 2 2x\ + 5a;| + z 4 
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Lemma 48. Suppose R is a square matrix with order n, P is an invertible square matrix 
with order k < n, Q G R kx ( n ~ k ), M G R( n ~ k )* k and N G r("-*0x(«-*0. jj R mn be 


written as partitioned matrix 


R = 


P Q 
M N 


then 

det(F) = det(P) det(7V — MP~ 1 Q). 


Proof. It is a well known result in linear algebra. □ 


For a square matrix M, we use to denote the determinant of the sub-matrix 

obtained by deleting the z-th row and the j-th column of M. 

Theorem 49. Suppose f G Z[z s ][x n ] is defined as in (3). Set F(x n ) = f(xf,... , a; 2 ). If 

(1) F(x n )\ Xi =o is nonzero and squarefree for any i € {1,... ,n}; 

(2) det(Aj) = i s nonzero and squarefree for any sub-sequence I o/{l,..., n}, 

and gcd(Aj 1,1 \ ..., A^ 7 ^ 7 ^) = 1 for any any sub-sequence |/| >2 of {1,..., n}; 

(3) gcd(P] 1,1 \ ..., Pj' 7 ^ 7 ^) = 1 for any sub-sequence |/| >2 of {1,..., n}; 

(4) det (Pi) is nonzero and squarefree for any sub-sequence I of {1,... ,n}; 

(5) gcd(det(P/),det(A/)) = 1 for any sub-sequence I of {1,... ,n}, 
then Np(P, [a;„]) = det(A„) det(A n+ i). 


Proof. We prove the theorem by induction on n. 

If n = 1, F{x i) = ai t ixf + 2a 1 ^xf + a 2 , 2 - Then 

OF 

Res(P, ^—, xfi) = 256a 7 1 a 2 , 2 (ai j ia 2 ,2 - a 2 2 ) 2 . 

OX\ 

By conditions (1), (2), (4) and (5), opi ^ 0,a 2j2 ^ 0, 01 , 102,2 — a\ 2 ^ 0, and a 2 , 2 and 
01 . 102,2 — af 2 are two coprime squarefree polynomials. Thus, Np(F, [aq]) = 01 , 1 ( 01 , 102,2 — 
a 2 , 2 ) = det(Ai) det(A 2 ). 

Assume that the conclusion holds for any quadratic polynomials with k variables 
where 1 < k < n. When n = k , let I be a sub-sequence of {1,..., n} with \I\ = n — 1. 
Without loss of generality, we assume that I = [1,..., n— 1]. Set A/ = (ojj)i<ij< n , B = 
( a l• • • 5 a n— 1 , 71)5 C = ( a l,n+l5 • • * 5 a n—l,n-fl)i and D = a n n X n 2a n n +\X n A O n _)_i, n _(_i. 
Then, F(x n ) could be written as 

n— 1 

P(*^n) ^ ' rtijX.j^Xj A ^ ) (az,n A ^n,i^)X^X n A (o n , n £ n A 2o n , n _j_iX n A O n -|_i,n-(-l) 

l<i,l<n i= 1 

=(x 1 ,. . . , X n _ 1: 1 )Pl(x ll . . . , X n _ 1 , 1) , 


where 


Pi 


Aj (Bxl A C) T 
Bx 2 n + C D 


By assumption, F\ Xi= o is squarefree for i G I , and det(Aj) is nonzero and squarefree, 
gcd(A^ -1) ,..., A ^ 11,17 |} ) = 1, gcd(P / ( , 1 ’ 1) ,.. .,Pjl 1 1,17 l} ) = 1, gcd(det(P//),det(A//)) = 1 
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for any sub-sequence /' of / with \I'\ < n— 1. Thus, by induction hypothesis, Np(P(a; n ), [*/]) = 
det (A/) det (Pj). 

In the following, we compute det(P/). By assumption, det(A/) = Pj n,n ' 1 is nonzero 
and squarefree. According to Lemma 48, 

det(P/) = det(Aj)(£> - (. Bx\ + C)Aj\Bx 2 n + C) T ) 

= det(Aj) ((a„ >n - BAj 1 B T )x* + 2(a„, n+ i - BA~ 1 C T )x 2 n + a„ + i, n+1 - CAj 1 C T ) 
= det(A/)(Aa;^ + 2/iz 2 + v), 


where A = (a„, n - BA T 1 B T ), fi = a„, n+ i - BAj 1 C T , v = a„+i,n+i - CAj 1 C T . By 
Lemma 48 again, 


det (A/) A = det(A/)(a n>n — BAj 1 B T ) 


= det 


Aj B t 

B CL n .n. 


_ ,(n+l,n+l) 

~ A n+1 > 

det(Ai)v =det(A/)(a„ +lf „ + i - CAj 1 C T ) 
\Ai C T 

C CL n - f-l,n+l 


= det 


_ A”’”) 

~- t± n +1 ’ 


(4) 


(5) 


Thus, both det(A/)A and det {Aj)v are the determinants of some principal sub-matrices 
of A n+ i with order n, respectively. 

Let H = det(P/), according to Lemma 48, it is clear that 


BH 

Res (H, — — ,x n ) =256 det(Aj) 7 A 2 istfj 2 — Xu) 2 

OXn 


=256 det (Aj ) 7 A 2 u det 


A /i 
l/i u 


=256 det (Ai ) 7 X 2 u det 


=256 det(Ai) b X 2 u det 


®n,n ^n,n+1 

| &n+l,n &n+l,n+l 

I Aj B t C t 

B dn : n ^n,n+1 

U c a n+l,n? &n+l,n+l 

=256 det(Aj) 5 X 2 u det(A n+1 ) 2 
=256det(A / ) 2 (A^’" +1) ) 2 4"+i ) det(A„ +1 ) 2 
=256(A^’")) 2 det(A„) 2 A^ } det(A„ +1 ) 2 . 

Since det (A/) and ^ n + 1 >”+ 1 ) are nonzero, according to (6), we have 



( 6 ) 
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Np(F, [x n ],n) =sqrfree(A^" | :" ) ^"’ Jl) det(A n ) det(A n+ i)). 

Similarly, for 1 < i < n, we have 

Np(F, [® n ], i) =sqrfree(A^A^’ l) det(A„) det(A„ +1 )). (7) 

By assumption, gcd( , • • • ,A[ l " ) l” ) ) = 1, gcd(A^ 1 ’ 1) , • • • 1 A { n ’ n) ) = 1, and det(A„), 
det(A n+1 ) are two nonzero squarefree polynomials with gcd(det(A„), det(A„ +1 )) = 1, 
thus 


Np(F, [aj„]) = gcd(Np(F, [x n \, 1),... ,Np(F, [x n ],n)) = det(A n ) det(A n+1 ). 
That completes the proof. □ 


Theorem 50. If the coefficients aij of f{x n ) in Theorem 49 are pairwise different 
real parameters except that aij = aj^, and F(x n ) = f(x 2 ,... ,xf). Then all the four 
hypotheses (l)-(5) of Theorem 49 hold. As a result, Np(F, [x. n ]) = det(A„) det(A„ +1 ). 


Proof. It is clear that the hypotheses (1) and (2) of Theorem 49 hold. We claim that 
that for any given to, |/| = in, P and clet(Pj) are pairwise nonconstant different 
irreducible polynomials in Z[a,ij\[x n \ for all n > to, so the hypotheses (3),(4) and (5) of 
Theorem 49 also hold. Here we denote a tJ = (ay i,..., a„+ i, n +i). 

We only prove that det(Pi) is a nonconstant irreducible polynomial. The other state¬ 
ments of the claim can be proved similarly. 

We prove the claim by introduction on n. If n = to, it is clear that the claim is true. 
Assume that the theorem holds for integers to < l < n — 1. We now consider the case 
l = n. Without of loss of generality, we assume that I = [1, 2,..., to]. 

Recall that F(x n ) could be written as 

F(x n ) = (x 2 m , l)Pp F (x 2 m , 1) T , 


where 


Pi,f = 


AI (pi,n? ? Pm,n ) 

(Pl,ny * iPm,ri) Pm-\-l,n 


and p i>n = YJj=m+i a i,i x2 j + a i,n+i for 1 < i < TO, and p m +i,n = S"= m +i a j,j x j + 
X/j=m+l 2 a n+l,j x j + ~l2m+l<i<j<n ^ a i,j X i X j F a n+l,n+l- Let B = (ay n , . . . , a m} n), C = 
(Pl,n— 1; ■ • • iPm,n— l)? and D Pm+l,n — a n ^ n X n T 2a nt n + iX n T 2 ay n x^ X n + 

a' n+l n+1 , where a' n+l n+1 is a polynomial with deg(a’ n+l n+1 ,x n ) = 0. Now Pi F can be 
simply written as 


Pr.F 


A, (Bx 2 n + C) T 
Bx 2 n + C D 
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In the following, we compute detBy Lemma 48, 
det(iV) = det(A/)(.D - (Bx 2 n + C)Aj\Bx 2 n + C) T ) 

= det (A/) ((a n , n - BAj 1 # 7 )^-!- 2(a n>n+ i + ^ a;,n£? - BAf 1 C T )x 2 n 

m+l<i<n 

+«; +1 , n+1 - caj 1 ^) 

= det (A/) (+ 2gx^ + v ), 


where A = (a„ )Tl — BAj 1 B T ), /.i = a„ ; „ + 1 +Sm+l<i<n a i,n x i BAj C T , V ^n+1,n+1 
CA~f l C T . By Lemma 48, 


det (A/) A = det(A/)(a n>n — BAj 1 B T ) 


= det 


A/ B t 

B Or, f. 


det(Aj)i/ =det(A/)(o^ +lin+1 - CA, 1 C ,T ) 

U/ C T 


= det 


C 


n+l,n+l 


We have det(A/)A = det(P/ ; G), det (Aj)u = det(P/ i //), where 


G{xm) — (x m , 1) 


A/ P T 
B a nn 


(^mi 1) j 


and 


B(x n — i) — (x m , 1 ) 


A/ C T 

C a n+l,n+l 


( X rri ; 1) • 


By induction, det(A/)A and det(A/)i/ are two different non-constant irreducible polyno¬ 
mials. Since deg(det(A/)/x,a„ ; „ + i) > 0, and deg(det(A/)A, a n , n+1 ) = deg(det(A/)i/, a n>n+ i) 
0, it is clear that det (Aj)fi ^ ±(det(A/)A • det(A/)i/ + 1), det(A/)/z ±(det(A/)A + 
det (A/)*/). Now the result follows from Lemma 51. We are done. □ 


Lemma 51. Let 1Z be a UFD with units ±1. Let a,b,c € 1Z, where b y^= ±(ac + 1), 
b ytz ±(a + c), and a, c are two non-unit coprime irreducible elements in TZ, then T{x) = 
ax 4 + bx 2 + c is an irreducible polynomial in TZ[x\. 

Proof. Otherwise, we may assume T(x) = g(x)h(x), where g,h are two nonconstant 
polynomials in TZ[x\. Notice that if a € TZ is a root of T(x), then —a is also a root of 
T(x), thus ( x 2 — a 2 ) is a factor of T. Thus, we may assume that deg(g) = deg(/i) = 2. Let 
9 = 9o + gix + 92 X 2 , h = ho + h±x + h 2 X 2 , where gi , hi £ TZ. By comparing the coefficients 
of T with gh, we have c = goho , 0 = goh± + giho, hig 2 + h 2 gi = 0. Assume that c\go, 
then c\ ho. if hi y^ 0, let l be the largest integer such that c l \h±, then l + 1 is the largest 
integer such that c l+1 \g\. But h\g 2 + h. 2 gi = 0, so c\g 2 , and c\gh, which contradicts with 
(a, c) = 1. We must have h\ = 0, and g\ = 0. We assume that go = c, ho = 1. Now, 
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there are four cases (g 2 = ±a,/i 2 = =Fl) or (g 2 = ±1 ,/i 2 = if a). All the four cases will 
contradict with the assumption that b ^ ±(ac + 1), b ^ ±(a + c). □ 

Theorem 52. Suppose g(x n ) = X!i<i j< n a i,j x i x j = x n A n xf i s a quadratic polynomial 
where a^j are pairwise different real parameters except that a^j = Oj j ,(l < i,j < n). Let 
G{x„) = g(x\,.. -,x 2 n ), then Np(G, [®„]) = det(A„). 

Proof. By Theorem 50, we have 

Np(G, [x n ], n) = Bp(det(A„_i) det(A. n )x n ,x n ) = det(A„_i) det(A„). 

Therefore, Np(G, [x n ]) = gcd(Np(G, [x n \, 1),..., Np(G, [®„], n)) = det(A„). □ 

By similar method we can prove that 

Theorem 53. Suppose g(x. n ) = X^i<i j< n a i,j x i x j = x n A n xf is a quadratic form where 
a.i^j are pairwise different real parameters except that aij = 0^,(1 < i,j < n ). Let 
A n = Then Hp(g, [x n ]) = det(A„) = discrim(g, [x n ]), where discrim( 5 , [x n ]) 

is the discriminant of the quadratic form g, it is an irreducible polynomial in r L[a i j\. 

This theorem implies that, for a class of polynomial g , Hp(g, [a; n ]) may coincide with 
its discriminant. 

Theorem 49 and Theorem 50 show that, for a generic copositive problem, we can 
compute the projection set Np(F, n — 1) = {/, Np(/, [, x n ]),..., Np(/, [x n ,..., x 2 )]} directly. 
Based on the theorem, it is easy to design a complete algorithm for solving copositive 
problems. However, for an input f(x n ), checking whether f(x n ) satisfies the hypothesis 
(3) of Theorem 49 is expensive. Therefore we propose a special incomplete algorithm CMT 
for copositive matrix testing, which is formally described as Algorithm 7. 

Remark 54. In Algorithm 7, we do not check the hypotheses of Theorem 49. Thus 
the algorithm is incomplete. However, the algorithm still makes sense because almost 
all f(x n ) defined by Eq. (3) satisfy the hypotheses. On the other hand, for an input 
f{x n ), checking whether f(x n ) satisfies the hypothesis (3) of Theorem 49 is expensive 
but the other three hypotheses are easy to check. Furthermore, f{x n ) is degenerate 
when some hypotheses do not hold and such case can be easily handled. Therefore, when 
implementing Algorithm CMT, we take into account those possible improvements. The 
details are omitted here. 

Complexity analysis of Algorithm 7. We analyze the upper bound on the number 
of algebraic operations of Algorithm 7. 

We first estimate the complexity of computing det(Pfc) for 1 < k < n. Because the 
entries of the last row and the last column of Pk are polynomials with fc 2 terms and 
the other entries are integers, we expand det(Pfc) by minors along the last column and 
then expand the minors again along the last rows. Therefore, the complexity of computing 
det(Pfc), i.e., gk, is 0(k 2 (k—2) 3 +k 2 (n—k) 2 ). Since gk{xk, £*) is an even quartic univariate 
polynomial, the complexity of real root isolation for gifxk, ot) is 0( 1) and we only need 
to choose positive sample points when calling SPOne. That means SPOne(gk(xk, a), 1) 
returns at most 3 points. Thus the scale of O in line 13 is at most 3* _1 . The cost of 
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Algorithm 7. CMT 

Input: An even quartic squarefree polynomial F(x n ) £ Z[x n ], n > 1, with an ordering 
x n < %n-i ■ ■ ■ ~< x\ and a set Q of nonnegative polynomials. 

Output: Whether or not F(x n ) > 0 on M" 

1: if F £ Q then return true 
2: end if 

3: for i from 1 to n do 

4: if CWT(F(xi,X 2 , ■ ■ ■ ,x n )\ Xi —o,Q) =false then return false 

5: else Q ■■=QU(F(x 1 ,X 2 ,...,x n )\ Xi = 0 ) 

6: end if 

7: end for 

8: g n (x n ) := det(P„) > Recall that F(x n ) = {x\,. .. ,^_ 1; l)P n (a;f,... ,xl_ 1: 1) T . 

9: O := SP0ne(g„, 1) 

10: for i from 2 to n do 

11: S:=0 

12: g n —i-t-l(Xn—i+li‘‘‘>X n ) det(P n _^_|_^) 

13: for a in O do 

14: S := SU(afflSPOne(p„_ i+ i(a; n _j + i,a),l)) 

15: end for 

16: O := S 

17: end for 

18: if £ O such that F(a n ) < 0 then return false 

19: end if 

20: return true 

computing gk(xk,a) is 0(k 2 ) for each sample point a.. Then the cost of the “for loop” 
at lines 10-17 is bounded by 

n 

- 2) 3 + i 2 (n - if + i 2 ^" 1 )) = 0{n 2 3"). 

i—i 

In line 18 of Algorithm 7, the number of checking F(a n ) is at most 3". And the 
complexity of every check in line 18 is 0(n 2 ) since F has at most n 2 terms. Then, the 
complexity of line 18 is bounded by 0(n 2 3"). Therefore, the complexity of lines 8 — 19 
is bounded by 0(n 2 3”). 

The scale of the set Q is at most Yfk =o (k) = ^°’ cos ^ recurs i ve ca ll s is 

bounded by 

0(£ [lj3 k k 2 ) = 0(n 2 2 2n ). 

In conclusion, the complexity of Algorithm 7 is bounded by 0(n 2 2 2n ). 

Remark 55. By a more careful discussion, we may choose at most two sample points on 
every call SP0ne(gfc(:rfc, a), 1). That will lead to an upper bound complexity, 0(3”n 2 ). 

8. Examples 

The Algorithm HpTwo, Algorithm PSD-HpTwo, and Algorithm CMT have been imple¬ 
mented as three programs using Maple. In this section, we report the performance of the 
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three programs, respectively. All the timings in the tables are in seconds. 


Example 9. (Strzebohski, 2000) 


/ = ax 3 + (a + b + c)x 2 + (a 2 + b 2 + c 2 )x + a 3 + b 3 + c 3 — 1. 


Under the order a -< b -< c -< x, an open CAD defined by / has 132 sample points, while 
an open sample obtained by the algorithm HpTwo has 15 sample points. 


Example 10. (Han et al. 20141 

/ = x 4 - 2x 2 y 2 + 2x 2 z 2 + y A 


2y 2 z 2 + z 4 + 2x 2 + 2y 2 - 4z 2 - 4. 


Under the order z y >- x, an open CAD defined by / has 113 sample points, while an 
open sample obtained by the algorithm HpTwo has 87 sample points. 

Example 11. For 100 random polynomials f(x,y,z) with degree Figure 3 shows 
the numbers of real roots of Bp(/, [z,y]), Bp(/, [y,z]) and Hp(/, [y,z]), respectively. It is 
clear that the number of real roots of Hp(/, [y, z\) is always less than those of Bp(/, [z, y]) 
and Bp (/, [y,z]). 


- BP(BP(f,z),y)=0 

- BP(BP(f,y),z)=0 

- gcd(BP(BP(f,z),y),BP(BP(f,y),z))=0 



Fig. 3. The number of real roots. 


Example 12. In this example, we compare the performance of Algorithm HpTwo with 
open CAD on randomly generated polynomials. All the data in this example were ob¬ 
tained on a PC with Intel(R) Core(TM) i5 3.20GHz CPU, 8GB RAM, Windows 7 and 
Maple 17. 

In the following table, we list the average time of projection phase and lifting phase, 
and the average number of sample points on 30 random polynomials with 4 variables and 
degree 4 generated by randpolyC [x,y,z,w] ,degree=4)-l. 


2 Generated by randpoly ([x,y ,z] ,degree=8) in Maple 15. 
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Projection 

Lifting 

Sample points 

HpTwo 

0.13 

0.29 

262 

open CAD 

0.19 

3.11 

486 


If we get random polynomials with 5 variables and degree 3 by the command 
randpoly([seq(x[i], i = 1..5)], degree = 3), 

then the degrees of some variables are usually one. That makes the computation very 
easy for both HpTwo and open CAD. Therefore, we run the command randpoly([seq(x[i], 
i = 1..5)], degree = 3) + add(x[i] 2 , i = 1..5) — 1 ten times to generate 10 random poly¬ 
nomials with 5 variables and degree 3. The data on the 10 polynomials are listed in the 
following table. 



Projection 

Lifting 

Sample points 

HpTwo 

2.87 

3.51 

2894 

open CAD 

0.76 

12.01 

7802 


For many random polynomials with 4 variables and degree greater than 4 (or 5 variables 
and degree greater than 3), neither HpTwo nor open CAD can finish computation in 
reasonable time. 


A main application of the new projection operator Hp is testing semi-definiteness of 
polynomials. Now, we illustrate the performance of our implementation of Algorithm 
PSD-HpTwo and Algorithm CMT with several non-trivial examples. For more examples, 
please visit the homepag^jof the second author. 

We report the timings of the programs CMT, PSD-HpTwo, and DPS (Han et al. 20131, the 
function PartialCylindricalAlgebraicDecomposition (PCAD) in Maple 15, function Find- 
Instance (FI) in Mathematica 9, QEPCAD B (QEPCAD), the program RAGlitf^] and 
S0STQ0LS in MATLAB^jon these examples. 

QEPCAD and S0STQ0LS were performed on a PC with Intel(R) Core(TM) i5 3.20GHz 
CPU, 4GB RAM and ubuntu. The other computations were performed on a laptop with 
Inter Core(TM) i5-3317U 1.70GHz CPU, 4GB RAM, Windows 8 and Maple 15. 


Example 13. (Han 2011) Prove that 


F ( x n) = C}2x 2 i) 2 

i—1 


^^2 x i x i +1 > o, 

i =1 


where x n +i = x\. 

Hereafter “oo” means either the running time is over 4000 seconds or the software 
fails to get an answer. 


3 https://sites.google.com/site/jingjunhan/home/software 

4 RAGlib release 3.23 (Mar., 2015). The RAGlib has gone through significant improvements. Thus, we 
updated the timing using the most recent version. 

5 The MATLAB version is R2011b, SOSTOOLS’s version is 3.00 and SeDuMi’s version is 1.3. 
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n 

5 

8 

11 

17 

23 

CMT 

0.06 

0.48 

1.28 

4.87 

11.95 

PSD — HpTwo 

0.28 

0.95 

6.26 

29.53 

140.01 

RAGlib 

0.42 

0.76 

1.34 

3.95 

8.25 

DPS 

0.29 

OO 

OO 

OO 

OO 

FI 

0.10 

oo 

OO 

OO 

oo 

PCAD 

0.26 

oo 

OO 

OO 

oo 

QEPCAD 

0.10 

oo 

oo 

OO 

oo 

S0STQ0LS 

0.23 

1.38 

3.94 

247.56 

oo 


We then test the semi-definiteness of the polynomials (in fact, all G(x n ) are indefinite) 

G(x n ) = F(x n ) — iq 10 x n - 

The timings are reported in the following table. 

n CMT PSD-HpTwo RAGlib DPS FI PCAD QEPCAD 

20 1.81 3.828 0.59 oo oo oo oo 

30 5.59 13.594 2.01 oo oo oo oo 

Example 14. Prove that 

3m+2 3m+2 m 

B(x 3m+2 ) = ( Y x i) 2 - 2 Y X iJ2 X i+V+l - °> 

i— 1 i—1 j =1 


where X 3 m +2 +r — x r . If m = 1, it is equivalent to the case n = 5 of Example 13. This 


form was once studied in Parrilo (2000). 


m 

CMT 

PSD-HpTwo 

RAGlib 

DPS 

FI 

PCAD 

QEPCAD 

i 

0.03 

0.296 

0.42 

0.297 

0.1 

0.26 

0.104 

2 

0.56 

1.390 

0.36 

23.094 

OO 

oo 

OO 

3 

0.71 

9.672 

0.75 

OO 

OO 

OO 

OO 

4 

7.68 

OO 

0.87 

OO 

OO 

OO 

OO 


Remark 56. As showed by Example 12, according to our experiments, the application 
of HpTwo and PSD-HpTwo is limited at 3-4 variables and low degrees generally. It is not 
difficult to see that, if the input polynomial f(x n ) is symmetric, the new projection 
operator Hp cannot reduce the projection scale and the number of sample points. Thus, 
it is reasonable to conclude that the complexity of PSD-HpTwo is still doubly exponential. 


29 















Acknowledgements 

Thank M. Safey El Din who provided us several examples and communicated with us 
on the usage of RAGlib. 

The authors would like to convey their gratitude to all the four referees of our IS- 
SAC’2014 paper, who provided their valuable comments, advice and suggestion, which 
helped improve this paper greatly on not only the presentation but also the technical 
details. 


References 

Andersson, L.-E., Chang, G., Elfving, T., 1995. Criteria for copositive matrices using 
simplices and barycentric coordinates. Linear Algebra and Its Applications 220, 9-30. 

Basu, S., Pollack, R., Roy, M.-F., 1998. A new algorithm to find a point in every cell 
defined by a family of polynomials. In: Quantifier elimination and cylindrical algebraic 
decomposition. Springer, pp. 341-350. 

Brown, C. W., 2001. Improved projection for cylindrical algebraic decomposition. J. 
Syrnb. Comput. 32 (5), 447-465. 

Collins, G. E., 1975. Quantifier elimination for real closed fields by cylindrical algebraic 
decompostion. In: Automata Theory and Formal Languages 2nd GI Conference Kaiser¬ 
slautern, May 20-23, 1975. Springer, pp. 134-183. 

Collins, G. E., 1998. Quantifier elimination by cylindrical algebraic decomposition-twenty 
years of progress. In: Quantifier elimination and cylindrical algebraic decomposition. 
Springer, pp. 8-23. 

Collins, G. E., Hong, H., 1991. Partial cylindrical algebraic decomposition for quantifier 
elimination. J. Symb. Comput. 12 (3), 299-328. 

Hadeler, K., 1983. On copositive matrices. Linear Algebra and its Applications 49, 79-89. 

Han, J., 2011. An Introduction to the Proving of Elementary Inequalities. Vol. 123. 
Harbin Institute of Technology Press. 

Han, J., Dai, L., Xia, B., 2014. Constructing fewer open cells by gcd computation in cad 
projection. In: Proc. ISSAC’2014. ACM, pp. 240-247. 

Han, J., Jin, Z., Xia, B., 2013. Proving inequalities and solving global optimization prob¬ 
lems via simplified cad projection. arXiv:1205.1224. 

Hiriart-Urruty, J.-B., Seeger, A., 2010. A variational approach to copositive matrices. 
SIAM review 52 (4), 593-629. 

Hong, H., 1990. An improvement of the projection operator in cylindrical algebraic de¬ 
composition. In: Proc. ISSAC’1990. ACM, pp. 261-264. 

Hong, H., Safey El Din, M., Jul. 2012. Variant quantifier elimination. J. Symb. Comput. 
47 (7), 883-901. 

URL http ://dx.doi. org/ 10.1016/j.j sc .2011.05.014 

McCallum, S., 1988. An improved projection operation for cylindrical algebraic decom¬ 
position of three-dimensional space. J. Symb. Comput. 5 (1), 141-161. 

McCallum, S., 1998. An improved projection operation for cylindrical algebraic decompo¬ 
sition. In: Quantifier Elimination and Cylindrical Algebraic Decomposition. Springer, 
pp. 242-268. 

Murty, K. G., Kabadi, S. N., 1987. Some np-complete problems in quadratic and nonlinear 
programming. Mathematical programming 39 (2), 117 129. 


30 


Parrilo, P. A., 2000. Structured semidefinite programs and senrialgebraic geometry meth¬ 
ods in robustness and optimization. Ph.D. thesis, CIT. 

Safey El Din, M., Schost, E., 2003. Polar varieties and computation of one point in each 
connected component of a smooth real algebraic set. In: Proc. ISSAC’2003. ACM, pp. 
224-231. 

Strzebonski, A., 2000. Solving systems of strict polynomial inequalities. J. Symb. Corn- 
put. 29 (3), 471-480. 


31 



